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singlet fermions that have the Yukawa couplings, as well as some Majorana mass terms. The 
number of such fields and the scales of their Majorana masses are not known. Several independent 
observations point to the possibility that some of these singlets may have masses well below the 
electroweak scale. A sterile neutrino with mass of a few keV can account for cosmological dark 
matter. The same particle would be emitted anisotropically from a cooling neutron star born in 
a supernova explosion. This anisotropy can be large enough to explain the observed velocities of 
pulsars. A lighter sterile neutrino, with mass of the order of eV, is implied by the LSND results; it 
can have profound implications for cosmology. We review the physics of sterile neutrinos and the 
roles they may play in astrophysics and cosmology. 
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STERILE NEUTRINOS IN PARTICLE PHYSICS 



CO 

The name sterile neutrino was coined by Bruno Pontecorvo, who hypothesized the ex- 
istence of the right-handed neutrinos in a seminal paper [1], in which he also considered 
vacuum neutrino oscillations in the laboratory and in astrophysics, the lepton number 
Ph| violation, the neutrinoless double beta decay, some rare processes, such as /i — ► ey, 

and several other questions that have dominated the neutrino physics for the next four 
decades. Most models of the neutrino masses introduce sterile (or right-handed) neutri- 
nos to generate the masses of the ordinary neutrinos via the seesaw mechanism [2]. The 
seesaw lagrangian 

££ = J%M +N a (if N a - yaaHLaNa - -y N c a N a + h.c. , (1) 

where J2?sm is the lagrangian of the Standard Model, includes some number n of singlet 
neutrinos N a (a = 1, ...,«) with Yukawa couplings y aa . Here H is the Higgs doublet and 
La ((X = e,n,x) are the lepton doublets. Theoretical considerations do not constrain 
the number n of sterile neutrinos. In particular, there is no constraint based on the 
anomaly cancellation because the sterile fermions do not couple to the gauge fields. 
The experimental limits exist only for the larger mixing angles yO. To explain the 
neutrino masses inferred from the atmospheric and solar neutrino experiments, n = 2 
singlets are sufficient [0], but a greater number is required if the lagrangian ([B is to 
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explain the LSND J5|l, the r-process nucleosynthesis [6], the pulsar kicks 
matter |[loL[lll [l^. ll3|l . and the formation of supermassive black holes [14] 



The scale of the right-handed Majorana masses M a is unknown; it can be much greater 
than the electroweak scale [2], or it may be as low as a few eV [[5, 13, 15]. Even if some of 



the right-handed Majorana masses are much larger than others, for example, if some of 



theM a (a = 1, ...,«/) are smaller than 100 GeV, while some others (a = «/, are much 
greater than 100 GeV, both classes can have a non-negligible contribution to the active 
neutrino masses. Obviously, this does not contradict the usual decoupling theorems, 
because the heavy states decouple from all the physical processes at low energies, but 
they can still contribute to the values of the active neutrino masses if the corresponding 
Yukawa couplings are large enough. 

Are they natural? 

The seesaw mechanism [0] can explain the smallness of the neutrino masses in the 
presence of the Yukawa couplings of order one if the Majorana masses M a are much 
larger than the electroweak scale. Indeed, in this case the masses of the lightest neutrinos 
are suppressed by the ratios (H) /M a . 

However, the origin of the Yukawa couplings remains unknown, and there is no 
experimental evidence to suggest that these couplings must be of order 1. In fact, 
the Yukawa couplings of the charged leptons are much smaller than 1. For example, 
the Yukawa coupling of the electron is as small as 10 . One can ask whether some 
theoretical models are more likely to produce the numbers of order one or much smaller 
than one. The two possibilities are, in fact, realized in two types of theoretical models. 
If the Yukawa couplings arise as some topological intersection numbers in string theory, 
they are generally expected to be of order one jlq] . although very small couplings are 
also possible fl 1 VM . If the Yukawa couplings arise from the overlap of the wavefunctions 
of fermions located on different branes in extra dimensions, they can be exponentially 
suppressed and are expected to be very small fl 1 811 . 

In the absence of the fundamental theory, one may hope to gain some insight about 
the size of the Yukawa couplings using 't Hooft's naturalness criterion II 1911 . which 
states essentially that a number can be naturally small if setting it to zero increases 
the symmetry of the lagrangian. A small breaking of the symmetry is then associated 
with the small non-zero value of the parameter. This naturalness criterion has been 
applied to a variety of theories; it is, for example, one of the main arguments in favor of 
supersymmetry. (Setting the Higgs mass to zero does not increase the symmetry of the 
Standard Model. Supersymmetry relates the Higgs mass to the Higgsino mass, which is 
protected by the chiral symmetry. Therefore, the light Higgs boson, which is not natural 
in the Standard Model, becomes natural in theories with softly broken supersymmetry.) 
In view of 't Hooft's criterion, the small Majorana mass is natural because setting M a to 
zero increases the symmetry of the lagrangian ([I]) [|20l l5|l. 

One can ask whether cosmology can provide any clues as to whether the mass scale of 
sterile neutrinos should be above or below the electroweak scale. It is desirable to have 
a theory that could generate the matter-antimatter asymmetry of the universe. In both 
limits of large and small M a one can have a successful leptogenesis: in the case of the 
high-scale seesaw, the baryon asymmetry can be generated from the out-of-equilibrium 
decays of heavy neutrinos 112 ill , while in the case of the low-energy seesaw, the matter- 
antimatter asymmetry can be produced by the neutrino oscillations Q22Q. The Big-Bang 
nucleosynthesis (BBN) can provide a constraint on the number of light relativistic 



species in equilibrium |23L 24], but the sterile neutrinos with the small mixing angles 



may never be in equilibrium in the early universe, even at the highest temperatures fl 1 Oil . 
Indeed, the effective mixing angle of neutrinos at high temperature is suppressed due 
to the interactions with plasma |26|l . and, therefore, the sterile neutrinos may never 
thermalize. High-precision measurements of the primordial abundances may probe the 
existence of sterile neutrinos and the lepton asymmetry of the universe in the future p"Hl . 

While many seesaw models assume that the sterile neutrinos have very large masses, 
which makes them unobservable, it is worthwhile to consider light sterile neutrinos in 
view of the above arguments, and also because they can explain several experimental 
results. In particular, sterile neutrinos can account for cosmological dark matter dlOO, 
they can explain the observed velocities of pulsars 10, 0,0], the x-ray photons from their 
decays can affect the star formation rf28ll . Finally, sterile neutrinos can explain the LSND 
result [@,|30|,[3ll], which is currently being tested by the MiniBooNE experiment. 



EXPERIMENTAL STATUS 

Laboratory experiments are able to set limits or discover sterile neutrinos with a large 
enough mixing angle. Depending on the mass, they can be searched in different experi- 
ments. 

The light sterile neutrinos, with masses below 10 2 eV, can be discovered in one of the 
neutrino oscillations experiments [32]. In fact, LSND has reported a result |33|] . which, 
in combination with the other experiments . implies the existence of at least one sterile 
neutrino, more likely, two sterile neutrinos pi 130]. It is also possible that sterile neutrino 
decays, rather than oscillations, are the explanation of the LSND result Bill . 

In the eV to MeV mass range, the "kinks" in the spectra of beta-decay electrons 



can be used to set limits on sterile neutrinos mixed with the electron neutrinos [13411 



Neutrinoless double beta decays can probe the Majorana neutrino masses [35]. An 
interesting proposal is to search for sterile neutrinos in beta decays using a complete 
kinematic reconstruction of the final state Q36Q . 

For masses in the MeV-GeV range, peak searches in production of neutrinos provide 
the limits. The massive neutrinos v,-, if they exist, can be produced in meson decays, 
e.g. — > /i ± v,-, with probabilities that depend on the mixing in the charged current. 
The energy spectrum of muons in such decays should contain monochromatic lines 0~4ll 
at Ti = (mfj + — 2m K m^ - my.)/2m n . Also, for the MeV-GeV masses one can set 
a number of constraints based on the decays of the heavy neutrinos into the "visible" 
particles, which would be observable by various detectors. These limits are discussed in 
Ref. [H. 



STERILE NEUTRINOS IN ASTROPHYSICS AND COSMOLOGY 

Sterile neutrinos can be produced in the early universe, as well as in supernova explo- 
sions. The light sterile neutrino, consistent with the LSND result, is consistent with the 



existing bounds on the big-bang nucleosynthesis [23, 24] and large-scale structure, es 



pecially if the mixing lepton asymmetry of the universe is larger than the baryon asym- 
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FIGURE 1. Sterile neutrinos with masses 2-25 keV can explain the pulsar kicks if the mixing angles 
are large enough (as shown). In the region marked excluded (x-rays), the relic sterile neutrinos produced in 
neutrino oscillations via the Dodelson-Widrow (DW) mechanism would have a density inconsistent with 
the existing x-ray bounds. If the sterile neutrinos constitute all the dark matter, their masses and mixings 
should fall below the dashed line. Note that DW mechanism is not sufficient to produce enough dark matter 
for points on the dashed line: a large lepton asymmetry lfl2tl or a new production mechanism iRil l38ll is 
required. The Lyman-a bound for dark-matter sterile neutrinos produced at temperatures T > 100 GeV is 
m s > 2.7 keV (see text or Ref. [38] for discussion). The cosmological and the x-ray bounds do not apply 
if the universe was never reheated above T ~ MeV (331. 



metry Il25h . A heavier sterile neutrino, with mass in the keV range is an appealing dark- 
matter candidate, as discussed below. 



Dark matter in the form of sterile neutrinos 



The sterile neutrinos can be the cosmological dark matter 01 EL EL El. The 
interactions already present in the lagrangian ([TJ allow for the production of relic 
sterile neutrinos via the Dodelson-Widrow (DW) mechanism [fioll in the right amount 
to account for all dark matter, i.e. Q. s « 0.2, if one of the Majorana masses is of the 
order of a keV. We will denote the dark-matter sterile neutrino V s and its mass m s . By 



assumption, the mixing angles are small, and 



m s « Mi (2) 

sine . f> 

The mass and mixing angle are subject to the x-ray limits on the photons from the 
decays of the relic sterile neutrinos rf39ll . as well as the Lyman-a bound [|40ll discussed 
below (see Fig. [I]). 

As was mentioned above, the relic sterile neutrinos can decay into the lighter neutrinos 



and an the x-ray photons 114 111 , which can be detected by the x-ray telescopes [|39ll . 
The flux of x-rays depends on the sterile neutrino abundance. If all the dark matter is 
made up of sterile neutrinos, Q. s ~ 0.2, then the limit on the mass and the mixing angle 
is given by the dashed line in Fig. [Q However, the interactions in the lagrangian (OQ) 
cannot produce such an Q. s = 0.2 population of the sterile neutrinos for the masses and 
mixing angles along this dashed line, unless the universe has a relatively large lepton 
asymmetry I112ll . If the lepton asymmetry is small, the interactions in eq. © can produce 
the relic sterile neutrinos via the neutrino oscillations off-resonance at some sub-GeV 



temperature [10]. This mechanism provides the lowest possible abundance (except for 
the low-temperature cosmologies, in which the universe is never reheated above a few 
MeV after inflation H3]). The model-independent bound [38] based on this scenario is 
shown as a solid (purple) region in Fig. Q] It is based on the flux limit from Ref. ll3~9ll 
and the analytical fit to the numerical calculation of the sterile neutrino production by 
Abazajian Ii42l1 . This calculation may have some hadronic uncertainties 14311 . but they 
appear to be under control for the mass and mixing angle in the range of interest 

If the lepton asymmetry of the universe is relatively large, the resonant oscillations 



can produce the requisite amount of dark matter even for smaller mixing angles [|12|1 . 
for which the x-ray limits are weak. (The x-ray flux is proportional to the square of the 
mixing angle.) It is also possible that some additional interactions, not present in eq. (OQ) 



can be responsible for the production of dark-matter sterile neutrinos [44, 38]. We will 
discuss this possibility in more detail below. 

The x-ray photons from sterile neutrino decays in the early universe could have 
affected the star formation. Although these x-rays alone are not sufficient to reionize 
the universe, they can catalyze the production of molecular hydrogen and speed up the 



star formation (12811 . which, in turn, could cause the reionization. Molecular hydrogen is 
a very important cooling agent necessary for the collapse of primordial gas clouds that 
gave birth to the first stars. The fraction of molecular hydrogen must exceed a certain 
minimal value for the star formation to begin 14611 . The reaction H+H^H2 + y is very 
slow in comparison with the combination of reactions 

H+ + H -> H+ + 7, (4) 
H+ + H -> H 2 + H+, (5) 

which are possible if the hydrogen is ionized. Therefore, the ionization fraction deter- 
mines the rate of molecular hydrogen production. If dark matter is made up of sterile 
neutrinos, their decays produce a sufficient flux of photons to increase the ionization 



fraction by as much as two orders of magnitude 112 811 . This has a dramatic effect on the 
production of molecular hydrogen and the subsequent star formation. 

Decays of the relic sterile neutrinos during the dark ages could produce an observable 
signature in the 21 -cm background Q29Q. It can be detected and studied by such instru- 
ments as the Low Frequency Array (LOFAR), the 21 Centimeter Array (21CMA), the 
Mileura Wide-field Array (MWA) and the Square Kilometer Array (SKA). 



New physics at the electroweak scale, and the Lyman- a bounds 

One can ask whether the mass M ~ keV in equation ([I]) is a fundamental constant 
of nature, or whether it could arise from some symmetry breaking via the Higgs mech- 
anism. For example, let us consider the following modification of the lagrangian (OQ) 
following Ref. [Q: 

& = & +N a (ifdp)N a -y aa HL a N a -jSN c a N a + V(H,S)+h.c., (6) 

where Jzfo includes the gauge and kinetic terms of the Standard Model, H is the Higgs 
doublet, S is the real boson, which is SU(2)-singlet, L a (a = e,/i,T) are the lepton 
doublets, and N a (a = l,...,ri) are the additional singlet neutrinos. Let us consider the 
following scalar potential: 

V(H, S) = m\ \H\ Z + m 2 2 S 2 + X 3 S 3 + X HS \H\ 2 S 2 + X S S 4 + X H \H\ 4 . (7) 

After the symmetry breaking, the Higgs doublet and singlet fields each develop a 
VEV, (H) = vo = 247 GeV, (S) = vi, and the singlet neutrinos acquire the Majorana 
masses M a = h a v\. The mass of the S boson in after the symmetry breaking is m s ~ vi. 
The presence of the singlet in the Higgs sector can be tested at the LHC Il47n . 

This modification makes no difference in the low-energy theory, for example, in its 
application to the masses of active neutrinos. However, the new coupling opens a new 
channel for production of sterile neutrinos in the early universe. Indeed, if the couplings 
of S to H are large enough, while h < 10 , the S boson can be in equilibrium at 
temperatures above its mass, while the sterile neutrino with a small mixing angle can be 
out of equilibrium at all times. This is the case, as long as the annihilations NN — > NN, 
NN — > scalars, etc. are not fast enough to keep the sterile neutrinos in equilibrium. Now, 
since S is in thermal equilibrium at high temperatures, some amount of sterile neutrinos 
can be produced through decays S — > AW. The amount of sterile neutrinos produced this 
way is determined by the h coupling (and is independent of the active-sterile neutrino 
mixing angle): 



1.4 x 10~ s / V m 



s 



where £ is the change in the number density of sterile neutrinos relative to T 3 due to 
the dilution taking place as the universe cools. For example, in the Standard Model, the 
reduction in the number of effective degrees of freedom that occurs during the cooling 



from the temperature T ~ 100 GeV to a temperature below 1 MeV causes the entropy 
increase and the dilution of any species out of equilibrium by factor t, ~ 33. 
At the same time, the sterile neutrino mass is determined by the VEV of S: 

keV 

/i(S)~keV (S) 10 2 GeV (9) 

h 

Based on the required values of Q, s and the mass, we conclude that the Higgs singlet 
should have a VEV at the electroweak scale. In this case most of the dark matter is 
produced at temperature above 100 GeV. 

At a lower temperature, some sterile neutrinos are also produced via the Dodelson- 
Widrow mechanism. This mechanism cannot be turned off. However, for small mixing 
angles, left of the yellow line in Fig. 1, the dominant contribution comes from the S- 
boson decays. 

This has dramatic implications for the Lyman- a bounds because the relation between 
the mass and the average momentum is very different from the Dodelson-Widrow 
case. The Lyman- a forest Il40ll constrains the free-streaming length of the dark-matter 
particles, but the relation between this length and the particle mass depends on the 
production mechanism. One can approximately relate the free-streaming length to the 
mass m s and the average momentum of the sterile neutrino: 

^ 1Mvc (w)(Ae^) (10) 



15 T 



T=lkeV 



The sterile neutrinos produced in the 5-boson decays have an almost thermal spectrum 
at the time of production. However, as the universe cools down, the number of effective 
degrees of freedom decreases from g*(T s ) = 110.5 to g* (0.1 MeV) = 3.36. Then £ = 
g*(T s ) /g*(0. 1 MeV) « 33. This causes the redshifting of (p s ) by the factor t, 1//3 : 



110.5 



(Ps)(T<mev) = 0.767/ ' (11) 



1/3 



Comparing eq. (PTT)) with the DW case, one concludes that, as long as the population of 
sterile neutrinos is dominated by those produced at a high temperature (large enough h, 
small 6), the Lyman-oc limit changes from 10 keV to 

m, > 2.7keV (12) 

This lower bound is shown in Fig. 1. 



Sterile neutrinos and the supernova 

Sterile neutrinos with masses below several MeV can be produced in the supernova 
explosion; they can play an important role in the nucleosynthesis, as well as in generating 
the supernova asymmetries and the pulsar kicks. 



The keV masses of sterile neutrinos coinside with the Mikheev-Smirnov- 
Wolfenstein [48] (MSW) resonance in nuclear matter for typical momenta of the 
supernova neutrinos |7|]. The position of this resonance is affected by the magnetic field 
due to the D'Olivo-Nieves-Pal-Semikoz effect Q49Q, which plays an important role in 
generating the pulsar velocities, as discussed below. Since the sterile neutrinos interact 
with nuclear matter very weakly, they can be very efficient at transporting the heat in the 
cooling proto-neutron star, altering the dynamics of the supernova [501]. This could lead 
to an enhancement of the supernova explosion. An additional enhancement can come 
from the increase in convection in front of the neutron star propelled by the asymmetric 
emission of sterile neutrinos rf52|l . 

In addition to playing an important role in the primordial nucleosynthesis ||24l. I51Q . 
sterile neutrinos can affect the r-process and the synthesis of heavy elements in the 
supernova Q6Q . 



The pulsar kicks 

The observations of neutrinos from SN1987A constrain the amount of energy that 
the sterile neutrinos can take out of the supernova, but they are still consistent with the 
sterile neutrinos that carry away as much as a half of the total energy of the supernova. 
A more detailed analysis shows that the emission of sterile neutrinos from a cooling 
newly born neutron star is anisotropic due to the star's magnetic field 10, |8Q. The 
anisotropy of this emission can result in a recoil velocity of the neutron star as high 
as ~ 10 3 km/s. While both the active and the sterile neutrinos are produced with some 
anisotropy, the asymmetry in the amplitudes of active neutrinos is quickly washed out in 
multiple scatterings as these neutrinos diffuse out of the star in the approximate thermal 
equilibrium 115 311 . In contrast, the sterile neutrinos are emitted from the supernova with 
the asymmetry equal to their production asymmetry. Hence, they give the recoiling 
neutron star a momentum, large enough to explain the pulsar kicks for the neutrino 
emission anisotropy as small as a few per cent fljj |8fl. This mechanism can be the 
explanation of the observed pulsar velocities 02]. The range of masses and mixing angles 
required to explain the pulsar kicks is shown in Fig. [TJ 

The pulsar kick mechanism based on the sterile neutrino emission has several addi- 
tional predictions 

• the kick velocities are expected to correlate with the axis of rotation 

• the kick should last 10 to 15 seconds, while the protoneutron star is cooling by the 
emission of neutrinos, but the onset of the kick can be delayed by a few seconds, 
depending on the mass and mixing angles 080. 

• neutrino-driven kicks can deposit additional energy behind the supernova 



shock |p2L 15011 . an d they are expected to produce asymmetric jets with the 



stronger jet pointing in the same direction as the neutron star velocity 115211 



Some of these predictions can already be tested statistically using the pulsar data [54]. 



CONCLUSIONS 



The underlying physics responsible for the neutrino masses is likely to involve the addi- 
tional SU(2)- singlet fermions, or sterile neutrinos. The Majorana masses of these states 
can range from a few eV to values well above the electroweak scale. Theoretical argu- 
ments have been made in favor of both the high-scale and the low-scale seesaw mecha- 
nisms: the high-scale seesaw may be favored by the connection with the Grand Unified 
Theories, while the low-scale seesaw is favored by 't Hooft's naturalness criterion. Cos- 
mological considerations are consistent with a vast range of mass scales. The laboratory 
bounds do not provide significant constraints on the sterile neutrinos, unless they have 
a large mixing with the active neutrinos. The atmospheric and solar neutrino oscillation 
results cannot be reconciled with the LSND result, unless sterile neutrinos (or other new 
physics) exist. 

There are several indirect astrophysical hints in favor of sterile neutrinos at the keV 
scale. Such neutrinos can explain the observed velocities of pulsars, they can be the dark 
matter, and they can play a role in star formation and reionization of the universe. 

The preponderance of indirect astrophysical hints may be a precursor of a major 
discovery, although it may also be a coincidence. One can hope to discover the sterile 
neutrinos in the X-ray observations. The mass around 3 keV and the mixing angle 
sin 2 ~ 3 x 10~ 9 appear to be particularly interesting because the sterile neutrino 
with such parameters could simultaneously explain the pulsar kicks and dark matter 
(assuming these sterile neutrinos are produced at the electroweak scale). However, it is 
worthwhile to search for the signal from sterile dark matter in other parts of the allowed 
parameter space shown in Fig. 1. The existence of a much lighter sterile neutrino, with 
a much greater mixing angle can be established experimentally if MiniBooNE confirms 
the LSND result. 
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